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I. INTRODUCTION 



One of the most fascinating aspects of classical field theory is the existence of travelling 
localized solutions of the non-linear equations describing physical systems: these solitons find 
application in a very broad class of physical problems. From the quantum field theory point 
of view, (being exact non-perturbative solutions of the classical theory) they are believed 
to carry information about the non-perturbative structure of the quantum theory. Their 
particle-like properties have an intuitive and reasonable interpretation as bound states of 
the elementary excitations of the corresponding quantum field theory. For specific two- 
dimensional models this idea was confirmed by explicit computations in the framework of 
the semiclassical approximation |lj , of quantum inverse scattering || and of other systematic 
expansions while, in more recent years, conformal field theory techniques have led to exact 
results for the S— matrix [|]]. Solitons also play fundamental role in recent non-perturbative 
developments of QFT (for a recent review, see ||), since they are seen as the "duals" of the 
elementary excitations appearing in the original Lagrangian ||. 

Solitons made their appearance in a completely different context some years ago, namely 
in low-energy phenomenological applications to physical systems confined to a plane. An 
interesting class of gauge theoretical models, describing non-relativistic matter coupled to a 
Chern-Simons gauge field, was introduced by Jackiw and Pi H] to obtain a simple realiza- 
tion of non-relativistic interacting anyons. For an appropriate choice of the self-interaction 
potential, in the static case, one can reduce the Euler-Lagrange equations to the completely 
integrable Toda equation (in the non-Abelian situation) or Liouville (in the Abelian one) 
||, with well-known soliton solutions. The dynamics of the system corresponds, in this 
case, to a dimensional reduction in time: recently [|9|jll| there has been instead considered a 
reduction in one spatial dimension. The physical reason lies in the hope that some charac- 
teristic of the model, in particular its fractional statistics, can be maintained by the related 
one-dimensional excitations and that, by a suitable modification, chiral behaviour can be 
induced. These two features are in fact relevant in the phenomenological description of the 
edge states in the Fractional Quantum Hall Effect [ Unfortunately the theory obtained 
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in this way, originally proposed in ||, partially fails to achieve its goal, because as was 
shown in [ID| no statistical transmutation arises for non-relativistic matter. On the other 
hand, as first observed in |[L0|| , a novel and interesting soliton structure is present there, 
finding its origin in the gauge coupling and in the chiral modification. Our investigation is a 
direct follow-up of the work begun in [TO] , and it is particularly devoted to the classical and 
quantum analysis of the solitons appearing in the theory. Actually, after realizing that the 
symmetry property of the original model is shared by a larger class of 1 + 1 dimensional the- 
ories, we have chosen to study a slight generalization of the model of P,|TU] , that fits equally 
well into the dimensional reduction procedure. Surprisingly, as we shall see, quantum effects 
naturally lead one to include the more general interactions that we have introduced. As in 
the analogous 2 + 1-dimensional family, integrability appears only for a suitable choice of the 
initial parameters (as observed in 0]), relating our system to the Derivative Non-Linear 
Schrodinger equation (DNLS) [jT3[| . Nevertheless soliton solutions exist for any value of the 
coupling constants. 

We start in Sect. 2 by describing the dimensional reduction of the Jackiw-Pi model and 
its modification, based essentially on the introduction of a chiral boson. The gauge action 
turns out to be of B — F type and we show that the system is equivalent to a general family 
of non-linear Schrodinger equations that do not possess Galilean invariance. In absence of 
a self-interaction potential for the matter, the model has non-relativistic scale invariance, 
an important ingredient in understanding some features of the theory: a suitable potential 
respecting the scale invariance can also be added. We discuss the conserved charges and we 
derive, from symmetry considerations, some general properties of localized classical solution. 
In Sect. 3 we compute the solitons, first without and then including a potential, and we 
discuss their properties: in particular they appear to be chiral, existing only for a fixed sign 
of the total momentum, and presenting an interesting particle-like behaviour, inherited from 
the scale invariance. We also consider finite energy solutions with non-vanishing boundary 
condition at (spatial) infinity (dark solitons), that possess opposite chirality with respect 
to the previous ones. We show that conserved charges must be carefully defined when 
non-trivial boundary conditions are present and we derive the relevant energy-momentum 
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dispersion relation. In Sect. 4, we discuss the quantization of the theory, paying particular 
attention to the relation between the solitons and the quantum bound states. Since the 
system is non-integrable for a generic choice of the coupling constants, it is difficult to 
obtain exact results, and we have explicitly solved the Schrodinger equation only for the 
two-body problem. Nevertheless we have obtained a perturbative solution, consistent with 
the Bohr-Sommerfield quantization of the solitons (in the weak-coupling limit) and the 
exact expression of the trace anomaly of the theory, showing that classical scale invariance 
is destroyed, except at the fixed point of the renormalization group, represented by DNLS. 
Combining this result with the semiclassical quantization of the soliton, we conjecture the 
general form of the energy for the n-body bound-state. Finally , in Sect. 5 we present our 
conclusion and indicate future directions. 



II. NON-LINEAR DERIVATIVE SCHRODINGER EQUATIONS FROM A 
GAUGE THEORY AND THEIR CLASSICAL SYMMETRIES 

A nonrelativistic gauge field theory that leads to planar anyons |7| is the nonlinear 
Schrodinger equation, gauged by a Chern-Simons field and governed by the Lagrange density 

£ (2+ i) = ^re a ^A a F Pl + im*{d t + iA^ - |1 £ \{d { + zA)*f - V(**V) (1) 

Here \I/ is the Schrodinger quantum field, giving rise to charged bosonic particles after 
second quantization. A^ possesses no propagating degrees of freedom; it can be eliminated, 
leaving a statistical Aharonov-Bohm interaction between the particles. V(^*ty) describes 
possible nonlinear self- interactions, and can be a general polynomial in the density We 
notice that the above Lagrangian is invariant under Galilean transformations, due to the 
topological nature of the Chern-Simons action. When analyzing the lineal problem [|1(J], it is 
natural to consider a dimensional reduction of ([!]), by suppressing dependence on the second 
spatial coordinate, renaming A 2 as (mc/h 2 )B and redefining the gauge field as A x = A x 
and A = A + mc 2 /2hB 2 . In this way one is led to a B-F gauge theory coupled to a 
non-relativistic bosonic field in 1 + 1 dimensions: 
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Ai+D = ^Be^F„ v + iW\d t + iA )V - + iA x )t>\ 2 - V(V*V), (2) 

2k 2m 

where k = (h 2 /mc)K is dimensionless and we have neglected d x (B 3 /3Hk) since it is a total 

spatial derivative. Eliminating the B and fields decouples them completely, in the sense 

that the phase of \1/ may be adjusted so that the interactions of the \1/ field are solely 



determined by V, and particle statistics remain unaffected ||14|| . In fact the equation of 
motion obtained by varying B, F^ u = 0, implies that the gauge field A^ is a pure gauge 
and it can be reabsorbed in a phase redefinition of ^. In order to obtain a non-trivial 
theory even in absence of V it is quite natural to introduce a kinetic term for B, which, 
for example, could be taken to be in the Klein-Gordon form. In the following, however, 
we prefer a simpler expression that describes "chiral" Bose fields, propagating only in one 
direction. The reason why we commit ourselves to this specific choice is twofold: we hope 
both to induce a statistical transmutation and that the chiral dynamics is inherited by 
Having a field theoretical model that was both chiral and with anomalous statistics would be 



of great relevance in the description of the edge states ||11|| . As we shall see, while this model 
achieves the first goal (chirality), it fails the other one. On the other hand its structure is 
so rich that a detailed study of its properties has intrinsic interest. 



The Lagrangian density for the chiral boson is proportional to ±BB' — vB'B' [ITS 



(Dot/prime indicate differentiation with respect to time/space.) Here v is a velocity and the 
ensuing equations of motion for this kinetic term (without further interaction) are solved by 
B = B(x±vt) (with suitable boundary conditions at spatial infinity), describing propagation 
with velocity ^v. Note that BB' is not invariant against a Galileo transformation, which is 
a symmetry of £(i+i) and of B' B'\ performing a Galileo boost on BB' with velocity v gives 
rise to vB'B', effectively boosting the v parameter by v. Consequently one can drop the 
vB'B' contribution to the kinetic B Lagrangian, thereby selecting to work in a global "rest 
frame" . Boosting a solution in this rest frame then produces a solution to the theory with 
a B' B' term. In view of the above, we choose the 5-kinetic Lagrange density to be 

and the total Lagrange density is C = Lb + The ensuing equations of motion are 



h 2 

ih(d t + iA )V + —{d x + iA l f^> - V'V = 0, (4a) 

F 01 -A^ , (4b) 
B' - ft/ctf** = 0, (4c) 
B + HkJ = 0, (4d) 

with J = ^ 0&*(d x + iA^ - ty(d x - L4i)\I>*) . The integrabihty condition for the last two 
equation leads to the usual continuity equation 

dt(W) + dJ=0. (5) 

Now in terms of field = exp f i [ dy A\(y : t) +i [ dt' Ai(x°,t') — i—B(x° , £)] \I> eqs 

\ Jx° Jt° Kfl J 

(|a|) and (f4b| ) become respectively 

. h 2 « * 

iW t <p + dl<p - h\j<p - V'4> = 0, (6) 

2m 

Foi = \dS*&. (7) 

The latter gives the electromagnetic field as a function of the density </>*</> while the former 
encodes all the dynamical content of the system. In the following we shall be mainly inter- 
ested in the case in which the potential V((f>*(f>) is absent: this corresponds to a Schrodinger 
equation with a current (j) nonlinearity: 

j = -^(4>*dJ-M*). (8) 

This is to be contrasted with the familiar nonlinear Schrodinger equation, where the non- 
linearity involves the charge density 



mi = -^d 2 j-\(^*)l (9) 

Note that the gauge interaction has dynamically produced a non-trivial alternative to the 
usual non-linear Schrodinger equation, much studied both from physical and mathematical 
points of view. However eq. (§), m contrast to eq. @, does not possess a local Lagrangian 
formulation directly in terms of the field 0. A local Lagrangian can be instead exhibited for 
the equation 
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ihd t c 



h 



+ 



+ V'c 



(10) 



governing the gauge equivalent variable 
Consider, in fact, the action 



exp 



A ' rX dy(4>*4>)(y,t) -fdt'jix^t'] 



i— 



t° 



S 



dt dxC 



dt dx 



ih 



(fd t <P - <f>d t <j>*) 



2m 



V(p 2 



the Euler-Lagrange equations that follow can be easily shown to reproduce eq. fliTf). In eqs. 
(|l0l) and ([TT]) p 2 represents the density 0*0, while 

A 



(12) 



is the corresponding current. 

The invariance of the action 5* under space/time translations reflects itself into the pres- 
ence of a divergenceless energy momentum tensor T A : 



T° = H 



2m 



2m 



(13a) 



T° x = V = mJ — ft~p 4 



T% = - \d x <p\ 2 - ^ay + V'p 2 

m Am 



V. 



(13b) 



where D x stands for the "covariant" derivative d x + i\/2p 2 . For solutions obeying suit- 
able boundary conditions, we can therefore write a conserved momentum P = J dx V and 
energy E — J dx 7i. The unconventional form of the momentum density signals lack of 
Galilean invariance. In fact, in contrast to the Galilean invariant case, V is not propor- 
tional to the U(l) current J. Indeed computing the time derivative of the Galileo generator 
G = tP-mfdxx?,*e find 



dG f A f 

—j^ = J dx(V — mJ) = — h— J dxp A , 



(14) 



namely G, depending on the sign of the coupling constant, always increases or decreases in 
time. On the other hand, for a particular choice of the potential V, the theory becomes 
scale invariant. In fact the action (|IT|) is unchanged under a dilation, t — > a 2 t, x — > ax, and 
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4>(x,t) — > ai(J)(a 2 t, ax), iff V(p 2 /a)a 3 = V(p 2 ): V(p 2 ) oc p 6 . The generator D of the scale 
symmetry takes the form 



D 



dxV = tH — - I dxxV, 



(15) 



where the density D = tT° — \xT\ obeys the continuity equation 



d t v + d x (t n - \x n - ^d xP 2 



-0. (16) 
We can remove the last term in eq. (|T6|) proportional to the derivative of p 2 by adding a 



superpotential to the energy-momentum tensor. In fact if we define an improved T y 



1 — 1 o™ U xP ■> 



8m 



Ti 2 

rpx rpx _ __£)2 T rpG 



rpx rpa 



(17) 



we obtain 



d t V + d x [tT* o --xT x x )=0. 



(18) 



The new energy momentum tensor satisfies 2T° = T x x , which is the non relativistic analog 
of the relativistic traceless condition for scale symmetry. In a Galilean invariant theory this 
condition would also ensure conformal invariance; here, instead, the absence of Galilean sym- 
metry requires the absence of conformal symmetry. Notice that the classical symmetries of 
our system form a three dimensional Lie sub-algebra of the conformal group (non-relativistic 
Weyl group in (1 + 1)— dimensions) with Poisson brackets 



{H,P} = 0, {D,P} = P, {D,H} = 2H. 



(19) 



In the following we shall be interested in classical "localized" solutions of the scale in- 
variant Hamiltonian 

2 



h 2 . 

H = — I dx 

2m 



dx + i^p 2 



4 



*\3 



(20) 



namely normalizable field configurations (N = J dxp 2 < oo) with finite energy and momen- 
tum. Here £ is a dimensionless parameter governing the strength of the cubic potential. 
In the absence of V (i.e. £ = 0) we recover the model of ref. P,[l0| that will be our main 



concern and whose solutions we shall present in detail. Some results about £ 7^ will be 
also discussed as well because of their relevance in the quantum theory. 

Unlike the conventional nonlinear Schrodinger equation, the family of scale invariant 
equations we have described, and in particular the one without cubic potential, does not 
appear to be completely integrable and thus analytic expressions for multisoliton solutions 
are not available. It has been remarked, however, that for £ = 1 the resulting Schrodinger 
equation becomes an integrable nonlinear, derivative Schrodinger equation with nonlinearity 

i^pdj 0. 



III. CLASSICAL SOLUTIONS 

When it comes to the problem of solving non-linear differential equations, very few 
general tools are available. One possibility is to look for solutions that possess particular 
symmetries or whose specific functional dependence simplifies the structure of the original 
equation. In our case a simple Ansatz, which allows us to integrate eq. (|T0D, is to assume 
that the density is a function only of x — vt, i.e. p(x,t) = p{x — vt). Drawing from the 
experience of integrable models, we can expect that this choice will allow us to explore 
the presence of one-soliton solutions or more generally of traveling waves. Substituting the 
ansatz into the continuity equation, yields 

d x (-vp 2 (x -vt) + J(x,t)) = 0, (21) 

and hence 

J(x, t) = vp 2 (x - vt) + Joo (t) . (22) 

Here Joo(t) is an arbitrary function of time, whose physical meaning will become transparent 
later. [ It is already clear that such a function is related to the boundary conditions chosen 
at infinity. ] In the following we shall be concerned with solutions that approach the vacuum 
at spatial infinity. Since TC is positive definite in the absence of the potential V(p 2 ), the 
vacuum solutions is first constrained by the requirement 
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H = =>- Ii=(d x + i^<f> = 0, (23) 

which is solved by taking <ft of the form 

(f)(x, t) = p (t) exp (-iu(t) - i^p 2 (t)x\ . (24) 

Requiring (]24]) to satisfy the equation implies that po(£) an d w (^) must be constants. Thus 

'A 



0(x, t) vacuum — Po 6Xp 



-z ( '-p\x + # 



(25) 



Let us notice that the vacuum solution (|25|) is characterized by a vanishing current J(x,t). 
Thus requiring that a solution approach the vacuum unambiguously fixes the value of Joo(t) 
in eq. fl22"|) , namely 

J(x,t) =i;(p 2 (a;-^)-p2). (26) 

This brief analysis of the space of vacua suggests that we distinguish two possible cases: 
solutions that approach the trivial vacuum <p vacU um = (po = 0) at spatial infinity, and 
the solutions that, instead, go to the ground state that has constant density po ^ 0. The 
former are characterized by well-defined energy, momentum and U(l) charge. The latter 
have finite energy, but infinite momentum and U(l) charge. This different behavior is easily 
understood in terms of the properties of the vacuum. In fact, while in both cases the vacuum 
is is invariant under time translation, the space translation and U(l) transformation leave 
it invariant only if po = 0. 



A. Solutions around the trivial vacuum (po = 0) 

This class of solutions is strongly constrained by the symmetry of the problem. Let us 

show how a wide number of their properties can be derived without explicitly solving the 

equation of motion. To begin with, the simple Ansatz p = p(x — vt), by means of eq. (PB|). 

implies that the momentum density V is a function of x — vt as well. Thus the dilation 

charge takes the form 

If 00 v r°° v 1 

D = tH - - / dx{x - vt)V{x - vt) --t dxV{x - vt) = t{H - -P) - -D , (27) 

2i J — oo 2i Jog £ £ 
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/oo 
dxxV(x). Since D is conserved and consequently time-independent we 
-oo 

obtain 

H = V -P, (28) 

namely the dispersion relation of a non-relativistic particle. Note that the result holds 
independently of the value of £. In fact this reasoning does not make use of the explicit form 
of the Hamiltonian, but only of its symmetries. 

To study further properties, we introduce the "center of mass" coordinate 

/oo 
dx xp 2 (x, t) 

.'V.uin =-M • (29) 

dx p 2 (x, t) 



This name is easily understood if we think of p 2 as the mass density. Its velocity will be 



L dx J(x,t) 

vcm = xcuit) = — - — 77 , with N = dx p 2 {x,t). (30) 

IV J— OO 

Here we have used the continuity equation for the current to eliminate the time derivative 
of the density. In this language the violation of Galilean invariance (|T4|) takes a suggestive 
form, 

X(P - mNv GM (t)) = -h— dxp 4 {x,t)<0. (31) 

Z J — 00 

Being valid for all t, this implies 

XP < mNX mm{v OM (t)} (32) 

For ^ < (i.e. repulsive potential and thus energy positive definite), one can show an 

analogous inequality for the energy. In fact let us consider the following inequality 

2 



00 
dx 

-00 



+ W- — : D X (, 
Ami 



> 0, (33) 



where w is an arbitrary parameter. In terms of the physical quantities 



E 

N + wNv CM {t) + w 2 -^- > 0. (34) 

2m 
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with Eq = — / dx \D x cj)\ . The fact that the previous equation holds for all w entails 

2m J-oo 



E > (35) 

but, for C, < 0,the real energy H is greater than Eq and thus 

E > mNv CM^)_ (36) 

Let us apply these inequalities to our Ansatz. A simple computation shows that v cm (t) = v. 
Thus eq. ( PTj ) can be rewritten as 

- -g— J < 0, (37) 

where we used eq. (f28l). For £ < 0, because of (|36|) , it implies At> < 0, i.e. the soliton is 
"chiral". In other words, given the sign of the coupling constant A, the sign of the velocity 
is determined. Recalling that 

A f°° 

P = m Nv -h- p 4 (x, t)dx, (38) 

2 J—oo 

the chirality can be also written in the form XP < 0. In this weaker form it is actually true 
for £ > as well (see below). 

We return now to the problem of finding explicit solutions. Upon introducing the 
parametrization 

<f>(x, t) = p(x — vt) exp[i0(x, £)], (39) 

eq. (P^| ) becomes 

TfLV A TfLV A r%—vt 

d x 9(x,t) = — --p 2 (x-vt) 0(x,t) = —x + e (t)--J_ go dyp 2 (y), (40) 

Substituting ([|(J) into the equation of motion, we obtain 

2m a , x , s „ , s m 2 v 2 , , 2mv , „ , . , , . 

—9 (t)p(y) = p"(y) - - -^-V(y)- (41) 

Here y = x — vt and the prime denotes the derivative with respect to y. The consistency of 



eq. plD requires that #o(t) is a constant cjo- Thus our equation can be rewritten as 
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2m / mv 



2mv 



-Xp 3 = 0, 



(42) 



or equivalently 



(p 



2„,2 



A2 



m v 



Ah' 



-IP 



mv 4 



where —ti—J 



. The possible arbitrary constant integration in eq. 



(43) 
is fixed 



8ft ' V" U ' 2h , 

by imposing that p — > as x — > ±oo. The problem of finding normalizable solutions is thus 
reduced to computing the zero energy orbits for a particle moving in a effective quartic 
potential. It is well-known that non-trivial (i.e. not identically constant) zero-energy orbits 
exist if and only if 



Xv < 



and 



7 > 0. 



(44) 



The first condition is particularly intriguing: given the sign of the coupling constant, a 
solution can be found only for a given sign of v, namely the system is "chiral". Integrating 
eq. fl4"3|) in the allowed parameter region, we obtain 

m\v\ 



p(x-vt) = - 1 





mv 




/ 


7 sech 




Xh 



2h 



-<fy(x -x - vt) 



(45) 



The phase can be in turn computed with the help of eq. (ffT|) 

■y/l(x -x - vt) 



= -^—x + u t + t) — - — signf A) tanh 
n 4 



m\v\ 



2Ti 



(46) 



The soliton's dynamical parameters, such as the U(l) charge, the energy or the momentum, 
can be now evaluted. Setting the U(l) charge 



N 



-oo A 



(47) 



we find that the energy and the momentum of the solitonic solution are 



E 



-Mv z 



and P = Mv where M = mN(l 



-X 2 N 2 ). 



(48) 



2 12 
The soliton's characteristics are those of a non-relativistic particle of mass M, moving with 
velocity v and composed of iV "constituents". Notice that iV is a function of v and u : for 
a given N we can get solitons of arbitrary velocity simply by tuning the phase velocity. 
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As we have seen in the previous section many properties are shared by the whole class 
of scale invariant systems: it is therefore interesting to explore the one-soliton structure of 
our theory when a potential cubic in p 2 is present. We parametrized V(p 2 ) as 

h 2 X 2 i 



5m 



-P 



(49) 



with £ e (— oo, +00). For £ = 1 we recover the integrable model described by the non-linear 



derivative Schrodinger equation, studied in fll2 |. Taking into accounts vanishing boundary 
conditions for p we simply have 



, 2 m 2 v 2 2 mv 1 



(p'Y - -^ip 



0. 



After defining Z 

(Z 
(Z 



m 2 v 2 



l\2 



2 2 

m v 



2 2 

,. TTI V n 

>\2 _ Z 2 



Atf ,r h " r ' 4 
p 2 and eq. ( |50D becomes 

( 7 e + 4)-(Z-2sign(At;)) 2 ] £>0 
( 7 |e|-4) + (Z + 2sign(At;)) 2 l £ < 0. 



(50) 



(51) 



Let us first consider £ > 0: following the same arguments of £ = case, we have normalizable 
solution for 7 > that are 



P 



mv 



7 



7^ + 4 cosh 



mv 



^7(x -x - vt) 



2sign(Af ) 



(52) 



In this range the solitons are not chiral, solutions exist for both the signs of the velocity: 
the condition on ujq is exactly the same as when £ = 0. Moreover it is immediately realized 
that the limit £ — > reproduces the correct solution only for Xv < 0, the other sign leading 
to a singular function. For £ < the situation is rather different: a short analysis of the 
"fictitious" potential yields the condition 

4 



Xv < 



and 



< 7 < 



lei 



(53) 



The functional form of p 2 is nevertheless similar, 



mv 

~x~h 



7 



4 — 7£ cosh 



mv 



y/l(x -Xq- Vt) 



+ 2 



(54) 



14 



This family of solutions is chiral and reproduces the solitons of the original equation. The 
different features of the two ranges and the particularity of the point £ = are clearly 
displayed when the conserved quantities are computed: 

4 



N 

E : 



|A|V? 



arctan , 



V4 + 7 £ + 2 sign(Aw) 
\| v /4+^-2 sign(Aw)' 



mv 
2\X\ 



(LlH nx - 2 sign w 



for £ > and 



N 



E = 



mv 



log 



(NWt\ 



(55a) 
(55b) 



1 + J2- V^ll 



1-J2- 



2|A| 



JV| A| 



tanh 



for £ < 0. We observe that for £ > iV takes value in a closed interval: 

7T 



(56a) 
(56b) 



< iV< 



7T 



< iV< 



|A|V? 
2tt 



Aw < 0, 
Aw > 0. 



(57a) 
(57b) 



No restriction arises for £ < 0. The formulae for the energy and the number present a 
potential singularity at £ = 0; actually for sign(Av) < the limit can be taken, recovering the 
previous results and the expressions in the two ranges are connected by analytic continuation. 
It is worth noticing that the energy is a non-polynomial function of iV and A (except for 
£ = 0), but an effective mass can be introduced as well, in perfect analogy with the zero 
potential case. 

Finally we observe that for all these solutions (for any £) the sign of the total momentum 
P is opposite to that of the coupling constant A, in spite of the fact that for £ > we can 
have negative and positive velocity (the total energy has the correct sign in order to realize 
the described situation). Chirality (in this more general meaning, i.e., \P < 0) is therefore 
a property shared by these scale invariant systems and it is clearly related to the non-linear 
current interaction present in the relevant Schrodinger equation. At the quantum level we 
shall recover such a property in the two-body quantum bound state and we shall make some 
observations for the general case. 
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B. Solutions around the non-trivial vacuum 



We investigate traveling solutions <f>(x, t) that approach, at spatial infinity, a vacuum 
with density po 7^ 0. When interpreting the standard non-linear Schrodinger equation as a 
description of a gas of bosons interacting via a repulsive two-body potential, this boundary 
condition is often called the "condensate boundary condition", since a time-independent 
solution p = po is usually understood as the Bose condensate. Finite energy solutions 
with the above boundary are therefore seen as travelling bubbles in a constant background 
density. On the other hand, the "naive" momentum P is not finite, because of the presence 
of the "unusual" p 4 term in the density T° x . In the presence of non- vanishing boundary 
conditions a conserved charge involves contributions coming from the boundary as well. 
The correct definition of the total momentum requires the addition of a suitable boundary 
term, that makes P different iable (in the functional sense) and therefore compatible with 
the Hamiltonian structure of the system. The resulting momentum turns out to be finite. 
This procedure entails a breaking of the relation E = |P, derived from the scale invariance, 
a fact that is nevertheless expected because our "condensate boundary condition" explicitly 
breaks such a symmetry. 

We start by presenting the phase of cf) when the non trivial boundary condition is taken into 
account: 




(58) 



and obviously it approaches to the phase of the vacuum as x — > oo. Then using eq.(^) we 
obtain the equation for p(y) 




(59) 



or equivalently (by using the boundary condition) 




(60) 



The relevant solution approaching to p\ as y — > oo exists if and only if 



Xv > 



and 




mv 



(61) 
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The shape of our "dark soliton" is 



p 2 (x-vt) = p\ 



mv 



Xhp, 



sech 



mvXpQ 

h 



mv 
UXpl 



(x — vt) 



(62) 



The dark soliton is chiral (Xv > 0) and moving in the opposite direction with respect to the 

soliton found in the previous section for vanishing cubic potential. Moreover only a finite 

range of velocity is permitted as clearly displayed by eq. (|6T|) : at variance with the soliton 

case we see that the only free parameter not fixed by the boundary conditions is the velocity 

hXpl 



v (we keep pi fixed); its absolute value cannot exceed the critical bound v r 
vacuum solution is recovered as v approaches v max _. 



The 



m 



Let us discuss the conserved quantities: the number, as defined in the previous section, 
is obviously not finite. It is natural, in the spirit of a bubble interpretation, to call number 
the quantity 

r+oo 



N 



/+oo 
dx (pi - p 2 
-oo 



that evaluated for eq.(^) is 



-1 



1. 



(63) 



(64) 



N is a monotonically decreasing function of v for any choice of background density and 
coupling constant. Next the expression for the energy is 

2 



E 



v 



The naive momentum P 



Vmax. 

X 



mJ + h-^p 4 



(65) 



is infinite, as we have anticipated, because 



of the presence of the "non-galilean" term (we recall that J — > at infinity). However 
our Schrodinger equation is a Hamiltonian system with phase space consisting of pairs of 
functions {4>{x), <f)*(x)) with boundary condition (<f)(x)<j)*(x)) — > p\ and the Poisson bracket 



{AB} 



-ih 



+oo 



dx 



5A 5B 



5A SB 



5<f)(x) S(j)*(x) 8(f)* (x) 5(f)(x) 



(66) 



To compute the Poisson bracket of P with some other functional we have to compute the 

5P 5P 

functional derivatives - , , , and „ , . , ; : using the expression of P as function the canonical 



5(f)(x) 5(f)* (x) 



pair we have 
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/+oo 
dx(54>d x (f>* - 8<j>*d x <f>). (67) 
-oo 



The functional derivative can be computed only if the first term on the right does not 
appear, so we have to add a functional that cancels the boundary term: this functional is 
proportional to 9' and the new momentum can be written in the compact form 

P = j J^dx(<p*d x <p - <R0*) (l - A.) . (68) 

P is automatically finite and obviously the subtraction is equivalent to a suitable improve- 
ment in the energy-momentum tensor. We define the improved components 

n = n _M^ ln (^, n=I ; + M a ,ln(l) (69) 

r+oo 

automatically satisfying the conservation equation and giving P = dxT° x . We remark 

J — 00 

that even the integral of T x x is convergent. A further observation concerns the scale invari- 
ance of the theory in presence of the condensate boundary condition: the generator D (that 
is not defined using the naive expression for P) is not conserved, as we can expect 

— = ^ / dxd t In — 70 

dt 2 7-oo \(f>* J 

The energy-momentum relation is therefore modified as follows 

v » dD 

"=2 P+ dT < 71) 

The momentum P of the dark soliton is easily evaluated 
P = -sign(A)ftpo 



2arctg,/— + h/( — -1] " \l— ~ 



^max. * . 4 / ^ \ / f^max. -, \ F^mo r. 



(72) 



that has been checked to be consistent with eq. (|71~1) . We notice that P is a monotonically 

decreasing function of v, ranging from +00 to (in absolute value), while the energy has a 

maximum value E = h-^— ( — f or v _ max - j s interesting to obtain the dispersion 

4|A| V 2 J 4 6 F 

dE dE 
relation for these solutions, namely — -. In a particle-like solution we expect — - = v. In 

dP dP 



ref. [16] the same subtraction procedure was followed in a Galilean invariant theory for dark 



soliton solutions: in that case the particle-like dispersion relation was recovered. This is 



not true in our case: performing the variation of the energy respect 0* and and using the 
equation of motion we get 



and therefore 



SE = v5P - v—p 2 8N, 



dE _ dP h\ 2 dN 
~dv ~ V ~dv~ V ^l pQ ~dv' 



(73) 



(74) 



Because the theory is not Galilean invariant the number depends on v, as we have seen, once 
Pq is considered a fixed parameter: the explicit computation gives 



dE _ (Vmax. ~ v){Vrnax. ~ 4f) 

dP 8v 2 - 10vv max . - v^ ax ' 



(75) 



which does not resemble the particle-like behavior. The curve E(P) can be easily plotted 



O . 2 



O . 1 




P 



E 



FIG. 1. The x-axis is — sign(A) — -s-, the y-axis is ., . 

The total energy goes to zero for small and large momenta and in particular as P — > we 
have 



E~ 2 - ^ \P\\ 



(76) 



3 hmpl 

We have seen that the finite-energy condition is equivalent to the requirement that J(±oo) = 
0. We can also imagine physical situations in which the asymptotic current is not zero; in 
this case a subtraction procedure can derived from the differentiability requirement, leading 
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to a sensible definition of the total energy. The general structure of the dark solitons in this 
case will be reported elsewhere. 

We turn now, as in the previous subsection, to the problem of generalizing the above 
analysis when a cubic potential in p 2 is added. A finite energy configuration requires that 
asymptotically (f> approaches a static solution with constant modulus and zero energy: the 
phase of these waves must satisfy 

(0' + \pl? = *U (77) 
Therefore only for £ > can we obtain the static, zero-energy solution 

= ~pg (l =F * + ft) (78) 

that carries the asymptotic current J = ± xfipt- When considering dark solitons inter- 

2m v 



polating these solutions at ±oo, this asymptotics reflects the following choice for J^: 

JUt) = -vp 2 ±^fU- (79) 
The relevant equation for the phase is easily solved 



while for the modulus we get 



, y2 , v2 mv (me A' 1 A 2 2 2 



(p') 2 = (P 2 -Plf 



(81) 



From the analysis of the above equation we can derive the parameter range for which solitonic 
solution exist: one has to solve a system of inequalities and here we state only the result, 
deferring details and physical interpretation to a forthcoming paper focused on "condensate 
boundary condition" solutions. 

a. ) The existing solutions are chiral with sign(At> ) > 0; 

b. ) There is an intrinsic lowest velocity 



2 



4/ilA., . 

M > o—tp'o 

6 m 
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There is no solution with finite energy unless < £ < 1; 

In this range we have two different kinds of solutions: setting 



4a 2 16a 2 16a 2 po 

Aa 2 Aa 2 ' 

16a 2 ' 
A = B 2 - AAC, 



with a = and /? = we find a true dark soliton 

2m 2m 



(82) 



2^4 

P 2 = Po(l~ 7=r rj. 77-^), (83) 

V A cosh(2p v — vtj) — B 

that represents an "hole" in the background density and a "bright" soliton 

p 2 =p 2 (l + ^ 2 ) (84) 

0V v / Acosh(2p oV / A(x -vt)) + B J V ' 

that is a positive density excitation over the background. Obviously while the first one goes 
smoothly to the solution for £ = the second one has a singular limit. The upper limit on 
the velocity is determined by the positivity of A. 

We end this section by remarking that "condensate boundary condition" solutions have 
no clear interpretation in the quantum theory; although their quantum role is an open 
problem they find applications in such branches of physics as condensed matter and plasma 
dynamics. 



IV. QUANTUM THEORY 

The simplest non relativistic example of the relation between solitons and quantum 
bound states is represented by the usual nonlinear Schrodinger equation. One may view 
eq. (|9|) as a Heisenberg equation for the quantum field <f)(x,t), which is taken to satisfy the 
commutation relation 
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[4>(x 1 ,t),4>*(x 2 ,t)] =S(xi- x 2 ). 



(85) 



The Hilbert space can be decomposed into invariant subspaces according to the integer 

/OO ^ a 

dx (f)*(f) and the n-body wave function, 
-OO 

<f> n (xi, -;X n ,t) = —j= < 0\<f)(xi,t)...(f)(x n ,t)\n >, (86) 



satisfies a Schrodinger equation with two-body, pairwise attractive ^-function interactions. 
The bound-state spectrum can be explicitly computed |18| and in the center of mass frame 
it is given by 

E n = -^(n 3 -n). (87) 
On the other hand, the semiclassical quantization of the one soliton solution at the leading 

dx 4 > *soliton ( f > soliton 

-OO 

with the quantum number n and the spectrum 

E S emici{n) = E C iass.(n) = -— r^™ 3 - (88) 

DAI 

Note that the classical energy of the soliton E C i aS s. gives the correct leading term as A — ► 
and Xn is kept fixed. Remarkably the next-to-leading correction diminishes n 3 by n (n 3 — > 
n 3 — n) and it is exact. From the quantum field theory point of view the correction linear in 
n to the energy stems from a renormalization effect, namely from the counterterm m(A)0*0 
necessary to cancel the divergences of the theory. 

One can argue on general ground @ and verify by computation ]l]) 12| , 17 1 that all inte- 



grable systems enjoy a similar stringent analogy between quantum bound states and solitons; 
much less is known for nonintegrable onesQ Our model (£ = 0) appears to belong to this 
second category: indeed it fails to pass the Painleve test. Nevertheless it can be consid- 
ered as a perturbed version of the derivative nonlinear Schrodinger equation (£ = 1) [p0[] , for 



1 From the point of view of the form factor approach, the question was recently investigated in |T!| . 
where a a non-integrable perturbation of the Ising model in external magnetic field was considered. 
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which some evidences of quantum integrability |21] exist. Thus it seems reasonable to expect 
that the quantization of the corresponding one-soliton solution will provide, in this case, a 
good approximation of the quantum spectrum. However before committing ourselves with 
the semiclassical approach, we shall try to extract some information through the pattern 
described at the beginning of this section. 

We take for the quantum hamiltonian the normal ordered expression 



H = I dx: 

2m 



d x + i^p 2 



(89) 



and we posit the canonical commutation relations Q8li|). The normal ordering prescription has 
been adopted to define properly the quantum Hamiltonian. This corresponds to removing all 
the singular interactions proportional to 5(0) in the resulting n-body Schrodinger equations. 
The above procedure is not without price, entailing in fact the loss of positivity for the 
quantum energy. [ The classical Hamiltonian is positive definite for any value of the coupling 
constant A.] With this choice, the n-body Schrodinger equation becomes 



2m 



n n y2 n 

Y^d 2 Xi + 2i\^25(xi- Xj)d Xt - — S(xi - Xj)5(xi - x k ) 

i=l i<j i^zj^zk 



E„ 



(90) 



[Xi, - -X ri 



cxp 



In eq. (p0|) the quantum non-integrability manifests itself through the presence of a three- 
body interaction generated by the (</>*0) 3 part of the Hamiltonian. In the case of the DNLS 
equation (£ = 1) this term is, instead, cancelled by an analogous contribution coming from 
the potential V. By translational invariance we can always separate the center of mass 
motion by introducing the parametrization 

n p x 1 n 

Yl — ~ Xn(xi,..x n ), with ^2d Xi Xn{xi, -X n ) = 0, (91) 

8=1 U J 1=1 

namely \ n depends only on the relative coordinates. However, because of the lack Galilean 
invariance, the total momentum P will not decouple completely and will appear as a pa- 
rameter in the reduced Schrodinger equation for Xn- 

From ( |90"D we note that the claim in of statistical transmutation for the quantum 
excitation of this model is inexact. It is impossible to remove the 5 interaction by a phase 
redefinition and therefore no change in the statistical behaviour should be expected |l(J . 
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A. The two-body problem 



We begin our quantum investigations by considering the simplest sector contained in the 
Hilbert space of the field theory defined by the quantum Hamiltonian (|89|): the two-particle 
one. Using eq. (|9T|) and defining x = x\ — £2 we find that \2 satisfies 

1" m + to + J = (92) 

The presence of the total momentum P in the 5— function potential for the relative motion 
vividly demonstrates the absence of Galilean invariance. Provided that 

\P < (93) 

eq. (|92]) possesses a bound state solution with energy 

E 2 = HL ( 1 _ X) = (94) 
2 4m ^ 4 / 2M RS . 1 ; 

Since P/2m can be identified with the velocity v of the center of mass, we recognize the 
condition ( p3|) as the one (see eq. (0)) that guarantees the existence of the soliton ( fi5D at 
the classical level. Remarkably even the scale invariance is preserved: in fact the dispersion 
relation (0) is unchanged under dilatation. [Unfortunately this property will fail in the 
n-body problem.] On the other hand the energy E2 and therefore the effective mass of the 
bound state 

2m . . 

become negative for A 2 > 4. The loss of positivity originates from the normal ordering 
adopted in defining the quantum Hamiltonian. As a three-body interaction of the form 
(<j)*(j)) 3 can never contribute to the two body problem, the bound state (|94j ) is, actually, 
shared by all the family (^) of models defined in sect. 2 and parametrized by the coupling 
constant £. Moreover all the solitons in this class of models enjoy the property XP < in 
perfect analogy with eq. (|93|) . 

It is instructive to recover the same results in a field theory framework, namely by 
resumming the perturbative series defined in terms of Feynman graphs. This approach will 
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be most useful in the discussion of the n— body case. The quantum propagator for is easily 
derived 



D(x,t) 



(27T) 



dk du 



(96) 



(h is put equal 1 from now on) and the interaction vertices are 





i XP/m 



-9i X (1- £)/2m 



FIG. 2. Feynman rules: P is the total incoming momentum. 



At one-loop level the relevant graph is 




FIG. 3. One loop contribution to the four point vertex. P is the total incoming momentum, k 
the loop momentum. 



The absence of the t and u channel implied by the non-relativistic character of the theory 
and moreover the absence of a tadpole contribution, because of the normal-ordering, greatly 
simplifies the final result. Choosing P 2 > 4mE, we get the one-loop amplitude 

A 2 P 2 1 



Ai(P,£,A) 



2m VP 2 - AmE 



(97) 



2 This restriction is the equivalent of the condition 7 > found in the discussion of the classical 
soliton in subsect. 3 A 
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Here P = pi + P2 is total momentum and E = uj\ + UJ2 is the total energy. To obtain the 
full amplitude we have only to resum the geometric series 



= xx>xxx + 



because there is no contribution from higher loops or from the six-vertex. [Recall in fact 
that the number of particles must be conserved in each stage of our diagrams since we are 
in a non-relativistic theory.] Thus we obtain 



A(P, E,X) = i 



,XP VP 2 - 4m£ 



m XP/2 + VP 2 ~4:mE' 



(98) 



p2 / 

We have a pole for XP < at the energy E = — 1 , recovering therefore the results 

4m y 4 J 

of the two-body Schrodinger equation. We remark that no infinity arises in the perturbative 
evaluation of the two-body scattering amplitude (once the normal-ordering is adopted) and 
no renormalization scale is needed: scale-invariance is maintained as the classical energy- 
momentum relation clearly shows. 



B. The n-body problem 



For the general n-body problem we did not succeed in computing the energy eigenvalues 
for the bound states exactly because of the presence of the three-body interaction that com- 
plicates the eigenvalue equation. Nevertheless starting from the usual non-linear Schrodinger 
equation we can construct a perturbative solution to the problem, that is reasonable for cer- 
tain values of the parameters and is in agreement with the soliton characteristics. The 
reduced equation is 



1 

2m 



n 2,XP n n 



1=1 



n 



i<3 



X 2 n 

— 5 ( X i ~ X j)$( X i ~ X k) 



Xn= \E n 



P 2 \ 



2mn 



Xr. 



Defining Z{ 



2XP 



Xj we obtain 



n 



(99) 



26 



]T -d\ + sign(AP) S(zi - Zj ) -A *sign(AP) £ - ^ + 



Wo 



A 2 n 

. • 

H 2 



Hi 



Xn 



ran 
2P 2 A 2 



P 2 
2mn 



Xn- 



(100) 



We recognize that the eigenvalue equation can be rewritten as the NLS one plus two per- 
turbations proportional to A and A 2 



A 2 

7io — A7Yi + ~X"^2 



(101) 



The NLS equation defined by Ho admits bound states only if AP < in agreement with the 
chiral nature of the soliton; thus we shall consider only this case. The zero-order problem 
-HoX^ = £<M } gives 



£f = -^(n 3 -n) 
y (o) = ,/ [2(n-l)]H 

and therefore the "zero-order" energy 



cxp 



"7 EN- 



(102a) 
(102b) 



2mn I 12 



PA 



(103) 



This result is correct as A — > 0, — is kept constant. Moreover the first correction, coming 

n 

from 7ii vanishes. Infact if £ n = £j® + A^ 1 ^ + .. the first order correction is simply 

/>oo n rco n 

e<p oc/ n? =1 ^(xi 0) )*(£<^ - %)^)xi 0) oc/ n« =1 ^( x ( ))* x (°)(Esign(^ - **)) = o, 

(104) 



because (Y^k s ^& n ( z i ~ z k)) — 0. Therefore 



E„ 



2mn 



1 - 



n 



Yin 



-A 2 +0(A 4 ). 



(105) 
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Obviously for n = 2 the result is exact: for generic n it shows that in a suitable limit the 

XP 

equation is essentially equivalent to NLS equation with effective coupling constant . It 

n 

is interesting to confront the perturbative solution with the exact energy spectrum of DNLS 



equation (7i 2 = 0), which has been recently computed |20| 

e dnls = ^_ A J*_ f _ t^l X 2 \ + 0(A 4 ). (106) 

4m / A\ 2mn \ I2n y ' y ' 

tan n arctan — 

V V 

The agreement between the perturbative and the exact result provides a strong check of our 
approach. 



C. Trace anomaly 

In subsec. 4A we have seen how to recover the two-body states by means of Feynman 
graphes. Here we would like to take advantage of the field theory approach to derive an- 
other exact result (connected to scale invariance) that would otherwise be quite difficult to 
establish. Classically the theory is invariant under dilatation; at the quantum level, instead, 
the ultraviolet divergences we encounter in perturbatively evaluating scattering amplitudes 
will spoil this property. In fact the well-known machinery of renormalization allows us to 
remove infinities consistently, but the price to be paid is the introduction of a "renormaliza- 
tion scale" that obviously breaks scale symmetry. In quantum mechanical language it is the 
highly singular three-body interaction that needs to be defined: different choices correspond 
to different self-adjoint extensions of a Schrodinger operator associate with the n— body 
problem. The (dimensional) parameter describing the possible extensions is usually related 
to the renormalized coupling constant ||22|| . 

A power counting analysis of the relevant diagrams shows that all divergences arise from 
the two loop six-point function or from its iterations (see figure below). No infinity comes 
from the 4-point vertex. 
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FIG. 4. Divergent 2-loop contribution to the six vertex function. P is the total incoming mo- 
mentum; k and 1 are the loop momenta. 



The two-loop amplitude in fig. 4. can be easily written down with the help of the Feynman 

rules given in sect. 4A. After integrating over the energy u% and 002 in the loops we are left 

( 9 

with the following integral over the momenta ( ao = A 2 (l — £o 



2m 

ial r +0 ° „ „ 1 



2m 2m 2m 

= -^2- / + " dk a i . (io7) 

Here, as before, E and P stand for the total energy and momentum respectively. [Although 
we are mainly interested in £ = case, we are forced to introduce a bare coupling constant 
£o 7^ in order to have a renormalizable theory] The final integral in eq. ( |107| ) is logarith- 
mically divergent. Upon introducing a cutoff A and a scale //, we can identify the singular 
contribution 

Af'{E,P) = ^-)og^. (108) 
12V37T fi z 

The divergent term can be reabsorbed in a redefinition of the coupling constant £ by defining 
a as follows 



eo = ^-^A 2 (l-^) 2 log^. (109) 



Since all the divergencies come from iterating the elementary graph in fig. 4 ( recall, in fact, 
that insertions of the four-point vertex do not give additional contributions to the divergent 
part), the result can be therefore made exact 
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I - 4n = rj- (HO) 



, x. A 2 



The /3-function obtained from this relation is 

= = -^ A2 (! - ^) 2 - ( m ) 

We remark that this result is exact and some of its consequences can be discussed. Firstly 
the coefficient of the three-body interaction appears to be running (it depends on the scale 
p 2 ): in field theoretical language it means that at the quantum level a potential V(p 2 ) ~ p 6 is 
induced from loop corrections. We are therefore naturally led to consider the general family 
parametrized by £. Secondly the fixed point of /?(£) is for £ R = 1, the DNLS equation, that 
is the integrable system: scale invariance is maintained in this model at quantum level as 
the exact expression, e.g., for the bound-state energy displays: 

„ P 2 „ , tanfiVarctan A) 

E = — , M = m i '-. 112 

2M' A v ; 

On the other hand for ^ ^ 1 we expect that the classical energy-momentum relation is 
broken at the quantum level. To make this idea more precise, we must relate the scale 
dependence of the quantum theory to the non conservation of the generator D of the dilata- 



tion. This can be done easily along the line of [23| where the same problem was discussed 
in 2+1 dimension for a p 4 theory: the methods relies in deriving a set of Ward identities 
from scale invariance. Following |23| we obtain for the proper vertices in momentum-space 
T n (pi, u>i, £r, //) the Ward identity 



77 

[d a + (3--)]T n (e- a Pi ,e- 2a cu i ,C R ,^ = -ir?(0, e~ a pi, e~ 2a uj h y), (113) 

T^(0, e~ a pi, e~ 2a u)i, £ R , /i) being the proper vertex derived from the Green function 

G^{y,x x ,..x n ) =< 0\T[(2T 00 (y) — T xx (y))(j)(x 1 )...(l)(x n )]\0 >, (114) 

while a is the dilatation parameter. If scale-invariance holds at the quantum level, the 
right-hand side of eq. ( |113j ) is zero. On the other hand the renormalization group equation, 
derived in the usual way is (we recall that no wave-function renormalization is present) 
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da d£ R 



n , 



(3-~)]T n (e- a p i ,e-"u i ,£ R ,ri = 0. 



(115) 



Comparing eq.( |113| ) and eq.( |115| ) we see that scale invariance is broken by the (3(£r) term. 
More precisely we can write the anomalous Ward identity 



d 



(116) 



that reflects in a quantum equation for the trace of the energy-momentum tensor 



The scale anomaly is 



rpa 



dD 
~dt 



2T 00 + 



327r m 



(117) 



V5* 



647T m 



+oo 



(118) 



From the above analysis it appears clear that physical quantities, like the bound-states 

d 

energies E n , must depend both on £i?(/i) and a so that u,—E n (£, fi) = 0. For example it 
is not difficult to compute the three-body bound-state in a simpler model, with potential 



purely proportional to p 6 {V(p 2 



(3 



° p 6 ; the family of models considered here reduces 



to this one in the limit A —>■ but A 2 £o = 2mao/9 is kept fixed). In this case, due to the 
absence of the two-body interaction, we have only to resum the geometric series 





leading to the exact scattering amplitude 

ia R 



A(P,E,a) 
that exhibits a pole for 



ni 



OtR 



12^ 



log 



P 2 - QmE 



P 2 p 2 Y2V3n 
E3= 6^-6^ eXph ^p) ] 



(119) 



(120) 



The second term, that is actually independent from /i 2 because the dependence of a R on /i 2 
is so chosen that this happens, is generated by the scale anomaly and it breaks the classical 
energy-momentum relation. In term of the physical parameters the scattering amplitude is 
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A(P,E,E 3 ) 



(121) 



m log 



P 2 - QE 3 m 
P 2 - 6Em 



that depends on one more dimensionful parameter (E 3 ). Moreover we see that a finite 
energy bound-state requires the bare coupling constant ao to be positive, implying that 
only the attractive potential leads to non-trivial physics, in close analogy with the analysis 
performed in p3| . In this particular case the non-trivial interaction can also be viewed as 
a self-adjoint extension of the two-dimensional laplacian, with one point removed. The self- 
adjoint extension parameter has the interpretation of the renormalized coupling constant 
a. 

In our case the resummation, even in the simplest case of the three-body scattering 
amplitude, is not a geometric series: it is not clear if only a bare parameter £o > 1 leads to 
a bound-state or if the derivative coupling (for XP < 0) can bind particles even in presence 
of a bare repulsive three-body potential. In the next subsection we will try to combine the 
soliton structure of the theory with the trace anomaly in order to get some information on 
the quantum bound-state of the theory. 



The semiclassical quantization of the one-soliton solution needs essentially the compu- 
tation of the small fluctuation around the soliton: following the approach of DHN JL], we 
have to calculate a functional determinant that produces the first quantum correction to 
the classical energy. In integrable systems, like NLS, DNLS or Sine-Gordon, one can use 
the multi-soliton solutions to compute the stability angles [|l|, from which the determinant 
can be determined. In particular, in all these systems, the contribution of the determinant 
consists in a (finite) renormalization of the parameters appearing in the initial lagrangian. 
In the DLNS case, e.g., the soliton has the mass Mci.{N) = while the mass of 



D. Solitons and bound states 



the quantum n— body bound-state is Mq {n) 



m tan[n arctan A] 



[p0|1 . The relevant renor 



arctan A 



malization is derived from the ratios 
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Mci.(N) taniVA Mq{n) tan[narctan A] , , 

/ — \~ = \ — ? and ~ — — = p ~r — , [iZZj 

Mci.(l) tan A MQ\n) tan|arctanAJ 

therefore consisting in the replacement A — > arctan A. [In general we must expect that only 
the mass ratios are correctly given by the semiclassical approximation.] In Sine-Gordon we 
have a similar renormalization for the coupling constant (3 2 

1 - fi 2 /8w 

The identification of the number parameter N, with the quantum number n, representing 
the number of particles bounded can be obtained applying the Bohr-Sommerfield rule 

dE = udn, (123) 

to the classical energy and taking uj the phase frequency at the maximum of the moving 
soliton (for x — vt — 0). 

In our case we do not have the two-soliton solution to use for the computation of the sta- 
bility angles and it is not clear what is the structure of the subleading quantum corrections. 
We can nevertheless compute the energy using the Bohr-Sommerfield formula, neglecting the 
contribution of the functional determinant. This can be easily done for the general family 
of scale invariant solitons, parametrized by £: we take the energies of the classical soliton 
(|55l) and ( ^6[) and we use the phase eq. (f40D , to obtain the frequency u that turns out to be 

mv 2 

lu = lu + —. (124) 

We rewrite the Bohr-Sommerfield formula as 

dE dn ,+ n „s 

7n=^7n' < 125 > 

that can be integrated to give n = N . The "zero" order energy for the n-body bound state 
is therefore 

E^(n) = E CL (n). 

Let us first study the case £ = 0; we expect that this equivalence is correct for A — >• 0. In 
fact 
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§§=i^ = ^-^ n2 - 1))+0{x% (126) 

Q- 1 -I 

12 

We stress that this relation is valid to 0(A 4 ) but in the case of n = 2 it is exact (the quantum 
corrections must compensate the higher terms in the A-expansion) . In the general case we 



see that eq. ( |126|) is in agreement with the perturbative expansion developed in Subsect.3. 
However we have seen that from renormalization the general scale-invariant potential is 
induced and we are naturally led to the ^-family of theories. We find that for £ > different 
answers are given, in the small- A limit, depending on the soliton chosen: 

= i(l + ^^(n 2 -l))+0(A 4 ), \v>0. (127) 

Only for £ = 1 do the two results coincide, while in general only for Xv < do we recover 
the "perturbative" computation: this is easily understood by realizing that only the soliton 
with Xv < is connected to an attractive NLS (that is the zero-order of the perturbative 
expansion). The meaning of the soliton with Xv > is not clear; moreover the entire 
procedure, except that for £ = 1, can be trusted only for P 2 — > oo and A — > 0, because 
we are neglecting the trace anomaly that drastically changes the momentum-energy relation 
respected by our "zero" -order approximation. Only in the large momentum limit does the 
classical term dominate over the "anomalous" contribution, as can be inferred from the 
analysis of the pure p 6 theory. 

In the following, in order to improve our knowledge about the quantum spectrum, we 
shall try to understand how the trace anomaly may modify the soliton analysis. We shall 
start assuming that the eigenvalue E n has the form 

E n = —F (x, 0?(/i 2 ), n, ^\ (128) 

The particular dependence on the mass m can be inferred from the n— body Schrodinger 
equation. The factor /i 2 has been pulled out to make the function F dimensionless. Since 
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the energy levels are physical observables, they must be independent of the subtraction point 
d 

//, i. e. /i— P n (/i) = 0, we have 



OF g(g OF 
dt 2 <9£ 



(129) 



P 



where t stands for — . It is not difficult to verify that the general solution for this equation 



is 



P 2 



A, n, — exp 



8tt 



,V3A 2 (l-£(/i)), 



(130) 



or equivalently 



P 2 



mM. 



A, n, — exp 



■ill 



,V3A 2 (i-e(/i)) 

The new dependence on the combination /i 2 exp ^ ^ X 2 {i-^{u)) 



;i3i) 



is just the symptom of the 



appearance of a new scale in the theory. This, for example, may be identified with the 
energy of the three-body bound state and it is obviously renormalization group invariant. 

In the limit P 2 — > oo, we expect a resurrection of scale invariance and consequently eq. 
( |131|) must match, in this limit, the soliton analysis. Taylor-expanding (|131|) we get 



E n 



P 2 



mMo [A, n] 



M.i[n, A]/i 2 exp 



57T 



,V3A 2 (i-e(/i)), 



O 



1 

P 2 



(132) 



The first term (the dominant one in the infinite momentum limit) respects scale invariance 
and we conjecture that the function 7Vfo[A,n] governing it is the same as in the DNLS 
equation. In fact we expect that such a model is the ultraviolet limit of our family of 
theories. An evidence of this is that £#(a0 flows to 1 for £r(/xo) > ^- Notice that in a 
Galilean invariant theory correction of order 1/P 2 are strictly forbidden, while they cannot 
be excluded in principle here. 

We conclude by observing that this brief analysis is consistent with the Galilean invariant 
model p 6 considered at the end of the previous section. 
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V. CONCLUSION 



In conclusion we have extensively studied a family of 1 + 1 dimensional theories that 
describes non-relativistic bosons interacting with a gauge potential: the gauge action has 
been chosen to be of B — F type plus a "crural" kinetic term for B. This form was suggested 
from the dimensional reduction of Chern-Simons theory coupled to non-relativistic matter 
and it represents a simple way to introduce chiral excitations that can be important in 
some condensed matter context. The theory was, in fact, proposed as relevant to modelling 
quantum Hall states: although it was shown in that it fails to achieve one of its goal 
(the statistical transmutation of bosons on a line), the system presents some very interesting 
features. First of all, it can be exactly reduced, solving for and B, to a self-interacting 
bosonic theory for which a local Lagrangian formulation is possible. Remarkably, the system 
possesses the non-relativistic scale invariance and for one choice of the parameters we recover 
an integrable equation (DNLS): for a generic choice we have a scale- invariant perturbation 
of the integrable model. The one-soliton structure of the theory has been examined, and 
it exhibits an interesting "chiral" behaviour: one-soliton solutions exist only for a fixed sign 
of the total momentum and they present a very peculiar particle-like energy-momentum 
relation inherited from the scale invariance. Solutions with non-trivial boundary conditions 
at infinity (dark solitons) were also found, existing for the opposite sign of momentum and 
with finite energy. In this case some care was needed in defining the conserved quanti- 
ties: in particular we have proposed a definition of the total momentum consistent with 
the Hamiltonian structure of the theory and leading to a finite result. The scale invariance 
is broken by our boundary conditions, which entails a complicate dispersion relation. Fi- 
nally we have studied the quantum dynamics, trying to discuss the relation between the 
classical one-soliton solutions and the quantum bound-states. The lack of integrability has 
not allowed a complete solution, but some result were obtained. The two-body problem 
was solved and the quantum bound states reproduce the chiral behaviour of the classical 
solutions. Moreover the perturbative computation of the energy for the N-body bound-state 
was found to be in agreement with the Bohr-Sommerfield quantization of the soliton in the 
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weak-coupling limit. On the other hands scale invariance was broken at quantum level by 
a trace anomaly: the fixed point of the renormalization group coincides with the integrable 
system, while, in the general case, we have deduced the functional form of the bound-states 
energy by combining the (classical) one-soliton solution with the informations coming from 
the (quantum) trace anomaly. There remains to prove or, at least, to check our proposal by 
explicit computations: for example to calculate the functional determinant representing, in 
the path-integral framework, the non-trivial quantum corrections or to solve explicitly the 
three-body Schrodinger equation. This subject together with the non-abelian extension of 
the model are currently under investigations. 
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